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The  purpose  of  this  study  is  to  evaluate  the  principle  of  the  technique 
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mine  the  propagation  direction  of  regional  seismic  phase  by  operating  on  a  com¬ 
prehensive  data  set. 

The  Smart  single-station  maximum-likelihood  surface-wave  processor  and  vari¬ 
ants  thereof,  all  relying  on  three-component  particle-motion  analysis,  demon¬ 
strated  in  this  study  their  utility  for  automatic  signal  azimuth  determination. 
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In  this  study  the  Smart  processor  was  found  to  be  that  which  simply  maxi¬ 
mized  the  ratio  of  energy  between  components.  Thus  the  P  wave  azimuth  is  deter¬ 
mined  by  minimizing  the  transverse  component,  the  Lg  by  maximizing  the  transverse 
and  the  emergence  cycle  by  finding  the  best  rectilinear  motion  fit  to  3  compo¬ 
nents  of  data.  For  the  entire  test  data  set,  the  average  azimuthal  error  was 
6.7  degrees  for  the  signals,  and  7.0  degrees  for  the  P  arrivals.  Combining 
the  estimated  L  azimuth  with  that  of  the  P  wave,  simply  by  taking  their  mean, 
increased  the  accuracy.  The  average  difference  between  the  mean  estimated  azi¬ 
muth  and  the  true  geodesic  azimuth  was  4.9  degrees  (7.0  degrees  rms) .  Moreover, 
the  F-statistic  computed  by  these  processors  serves  to  separate  poor  azimuthal 
estimates  from  the  population:  the  azimuthal  estimates  (over  80%  for  this  data 
set)  which  passed  the  arbitrary  F  threshold  set  for  this  study  differed  from  the 
known  geodesic  azimuths  by  an  average  of  3.9  degrees  (4.9  degrees  rms).  It  is 
the  author's  experience  that  such  accuracy  is  comparable  to  or  better  than  that 
from  a  well-desxgned  array  of  sensors.  .  i 

Besides  this,  the  results  of  the  study  suggest  that  particle-motion  analysis 
can  also  be  used  to  pick  seismic  phases  automatically  and  that  it  can  be  made  to 
yield  automatic  distance  determinations  simultaneously  with  the  azimuth  measure¬ 
ments.  One  complete  source-location  determination  could  thus  be  made  for  each 
station  recording  an  event. 
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ABSTRACT 


The  purpose  of  this  study  is  to  evaluate  the  principle  of  the  technique 
known  as  the  Smart  processor,  developed  at  the  SDAC  (Smart,  1977)  and  to 
determine  the  propagation  direction  of  regional  seismic  phases  by  operating 
on  a  comprehensive  data  set. 

The  Smart  single- station  maximum- likelihood  surface-wave  processor  and 
variants  thereof,  all  relying  on  three-component  particle-motion  analysis, 
demonstrated  in  this  study  their  utility  for  automatic  signal  azimuth 
determination. 

In  this  study  the  Smart  processor  was  found  to  be  that  which  simply 

maximized  the  ratio  of  energy  between  components.  Thus  the  P  wave  azimuth 

is  determined  by  minimizing  the  transverse  component,  the  L  by  maximizing 

S 

the  transverse,  and  the  emergence  cycle  by  finding  the  best  rectilinear 
motion  fit  to  3  components  of  data.  For  the  entire  test  data  set,  the 
average  azimuthal  error  was  6.7  degrees  for  the  signals,  and  7.0  degrees 
for  the  P  arrivals.  Combining  the  estimated  L  azimuth  with  that  of  the  P 

s 

wave,  simply  by  taking  their  mean,  increased  the  accuracy.  The  average  dif¬ 
ference  between  the  mean  estimated  azimuth  and  the  true  geodesic  azimuth  was 
4.9  degrees  (7.0  degrees  rms) .  Moreover,  the  F-statistic  computed  by  these 
processors  serves  to  separate  poor  azimuthal  estimates  from  the  population: 
the  azimuthal  estimates  (over  80%  for  this  data  set)  which  passed  the  arbi¬ 
trary  F  threshold  set  for  this  study  differed  from  the  known  geodesic  azimuths 
by  an  average  of  3.9  degrees  (4.9  degrees  rms).  It  is  the  author's  experience 
that  such  accuracy  is  comparable  to  or  better  than  that  from  a  well-designed 
array  of  sensors. 

Besides  this,  the  results  of  the  study  suggest  that  particle-motion 
analysis  can  also  be  used  to  pick  seismic  phases  automatically  and  that  it 
can  be  made  to  yield  automatic  distance  determinations  simultaneously  with 
the  azimuth  measurements.  One  complete  source-location  determination  could 
thus  be  made  for  each  station  recording  an  event. 
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INTRODUCTION 


This  study  was  prompted  by  the  modest  success  of  an  experiment 
(Smart,  1977)  in  which  a  single-station,  surface-wave,  particle  motion  pro¬ 
cessor  was  applied  to  certain  regional  seismic  surface-wave  data.  The 

signals,  combined  L  -R  waves  from  the  underground  nuclear  explosions  GNOME 
g  g 

and  SALMON,  were  recorded  at  several  Long  Range  Seismic  Measurement  (LRSM) 
stations  in  1961  and  1964  respectively.  The  processor,  which  fits  a  combined, 
fundamental-mode  Rayleigh  wave  and  Love  wave  to  given  3-component  seismic 
records,  would  seem  ill-suited  for  application  to  short-period  surface  waves 
recorded  at  distances  of  no  more  than  30  degrees  from  the  source.  However, 
in  spite  of  the  evident  mixture  of  higher  modes  in  the  records,  fairly 
reasonable  azimuths  were  estimated  with  the  processor,  although  it  could  not 
be  said  to  have  detected  the  signals. 

In  this  study  the  surface-wave  processor  is  applied  again  to  regional 
surface  waves,  and  another  processor,  suitable  for  first-arrival  P  waves, 
is  developed  and  applied  to  a  test  data  set  of  seismic  signals  recorded  at 
regional  distances.  The  objective  is  to  test  the  usefulness  of  such  single¬ 
station  particle-motion  processor  for  detecting,  identifying,  locating, 
and  analyzing  seismic  regional  phases. 

P-Wave  Processing 

Since  P-wave  particle  motion  is  often  quite  irregular,  departing  widely 
from  the  simple  motion  that  an  homogeneous  earth  would  produce,  (see 
Figure  1),  it  is  desirable  that  a  P-wave  processor  not  be  restricted  to 
fitting  data  to  a  model  of  straight-line,  back-and-forth  motion.  Here 
we  develop  a  process  which  merely  searches  out  the  orientation,  azimuth  and 
incidence  angle  in  which  the  mean  square  excursion  of  the  data  is  greatest, 
at  least  in  the  frequency  band  of  interest,  without  explicitly  modeling 
linear  motion.  However,  in  spite  of  our  maneuver,  this  development  is 
identical  to  fitting  a  straight-line  motion  model  to  the  data,  though  we 
do  not  explicitly  invoke  the  model  in  our  development.  This  identity  will 
be  demonstrated  after  the  development  of  the  processor,  which  follows  here. 

It  might  be  thought  that  this  processor  is  similar  to  "remode"  a  non¬ 
linear  processor  discovered  by  Flinn  in  Archambeau  et.  al.  (1965).  The 
present  processor,  however,  is  linear  and  from  this  difference  come  several 
features  which  make  it  superior  to  remode. 
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Development  of  a  P-Wave  Processor 
From  Figure  2 


r  =  zcos<|>  +  ysin4>cos0  +  xsin<{>sinO 


or 


r  =  yz  +  a  <5y  4-  g6x 


in  which 


and 


y  =  cos  <p 


6  =  sin  <p 


a  =  cos  0 


$  =  sin  0 


The  component  of  displacement  in  the  r  direction  at  the 
instant  is,  then 


.tk 

l 


sampling 


r .  =  yz .  +  a&y .  +  gdx, 

l  i  •'i  i 


The  sum  of  the  squared  excursion  along  the  r  direction  is 


N 

V  5  l 
i=l 


over  the  time  interval  i=l,...,N. 

To  move  into  the  frequency  domain  observe  that 


N 

v  =  l  r2 
i=l 


N 


N 

in  which  the  4 '  .  are  complex  coefficients  of  the  —  +  1  terms  of  the  discrete 
Fourier  transform. 
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Let 


in  which 


and 


and 

So 


and 


V  -  =  R,  +  iR' 

i  J  J 


Rj  =  real  [».] 

Rj  =  imaginary  [Yy] 


R .  =  yZ.  +  «6Y!  +  66X. 

j  J  J  J 

R'.  =  yZ'  +  a6Y’  +  B6X! 

J  j  J  J 


I  ¥  .I2  =  R2  +  R'.2 
1  J  J 


N 

N  2 

V  =  .  E,r2  =  E  [yZ.  +  a6Y.  +  36X . ] 2 
1=1  1  j=0  J  J  J 

+  [yZ^  +aBY^  +B6X^] 


Should  we  choose  to  do  so,  we  are  able  to  compute  the  mean  squared 
excursion  in  a  limited  band 


n 

V  =  E  R?  +  R'. 2 
m,n  .  j  i 
j=m  J 

in  which  m  and  n  are  arbitrary  integers  such  that 

0  <_  m  <  n 

m  <  n 

Then 

V  =  0262A  +  ct262B  +  YZC 
m,n 

+  a862D  +  By^F  +  »y5G 
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where 


A  = 

B  = 

C  = 

D  = 

F  = 

G  = 


E  X2  +  X ' 2 

j=m  j  j 


£  Y2  +  Y '  2 

j=m  j  J 


,£  Z;  +  Z 


i  2 


J=m  J  j 


2  E  X.Y.  +  X'.Y! 
j=m  J  J  J  J 


2  .  £  X.Z,  +  X’.Z! 
J=m  j  j  j  j 


2,£  Y  . Z  .  +  Y'.Z' 
j=m  J  j  J  j 


Dropping  the  subscription  on  v 

v  =  Y2b1  +  <Syb2  +  62b3 

in  which 

C 

6F  +  aC 

a2B  +  82A  +  agD 


Then 


2v  =  b^  +  b^  +  p(b^  -  b^)  +  5 b2 
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in  which 


?  =  sin  2<J>  and  y  =  cos  2$ 

Maximizing  v  with  respect  to  <f>,  the  incidence  angle,  for  arbitrary  azimuth  0, 

If  =~t(bi  -  b3)  +  ub2  =  0 

£  =  b2/den 


and 


in  which 


P  -  (b^  -  b2)/den 


den  =  /(b^  -  b^)2  +  b^ 


Maximized  with  respect  to  $,  then 


2  4-  K211/2 


2v  =  b  +  b3  +  [(b.^  -  b^j)^  +  b2] 


Now  maximizing  v  with  respect  to  azimuth,  0, 


2d v 
d0 


b5  +  I[(bl  "  b3)2  +  b21_1/2["2(bl  "  b3)b5  +  2b2b4]  =  ° 
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•* 


in  which 


b,  =  b'  =  aF  -  BG 
4  2 

and 

b5  =  =  2aB(A  -  B)  +  (a2  -  B2)D 

Then 

b5{(b1  -  b3)2  +  b|]1/2  =  (bL  -  b3)b5  -  b2b4 

b5(bl  -  b3)2  +  b5b2  =  (bl  -  b3)2b5  +  b2b4  '  2b2b4b5  (bl  ~  b3> 
b2(b4  '  b5>  =  2b4b5(bl  '  b3} 

This  can  be  written  in  the  form 

a3  mJ  +  a23x  +  cx32 1/  +  B3Z  =  0 

where 

W  =  (D2  -  F2)  G  -  2(B  -  C)DF 

X  =  4(A  -  B) (C  -  B) F  -  F[(F2-G2)  +  (D2-G2) ]  +  [ (A-B) 
y  =  4(A-B) (A-C)G  +  G[(F2-G2)  -  (D2-F2)]  -  2[(A-B)  - 
Z  =  (D2  -  G2)F  -  2 (A  -  C)DG 
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+  (A-C) ]DG 
(B-C) ]DF 


To  eliminate  g  from  the  equation,  note  that 


a3(w  -  y)  +  atj  =  -  Bla2(x  -  z)  +  z] 

and  square 

a6(w  -  y)2  +  y)  i/ +  a2t/2  =  g2(u2(x  -  z)  +  z]2 

=  a4(x  -  z)2  +  2a2(x  -  z)z  +  z2  -  a6(x  -  z)2 
-  2av(x  -  z) z  -  a2z2 

This  can  be  written  in  the  form  of  a  cubic  equation  in  a2: 

(a2)3P  +  (a2) 2Q  +  a2R  +  S  =  0 

where  now 

P-(w-.I /)2+(x-  z)2 
Q  =  2(m-  tji>(/+2(x-  z)  z  -  (x-  z)2 
R  =  </2  -  2(x  -  z)z  +  z2 
S  =  -z2 

The  solution  we  seek  is  restrained  by  these  considerations 

1)  0  <  a2  <  1  on  the  real  numbers 

2)  -  1  <  v  <  1 

and 

0  <  ?  <  1 


L 
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since 


,  TT 
0  £  4i  <  - 

from  the  constraint  that  the  particle  motion  in  the  positive  sense  must  be 
up  and  away  from  the  source  of  the  arriving  P-wave. 

The  cubic  equation  can,  of  course,  be  solved  analytically.  Thus  v  is 
optimized  without  a  computer  search,  which  ensures  processing  speed. 

To  show  that  this  development  is  identical  to  modeling  a  straight-line, 

back-and-forth  motion,  consider  such  motion  along  azimuth  6  and  incident 

th 

angle  <p,  having  waveform  r  where  r^  is  the  i  sample  in  the  time  series 
representation  of  the  model.  The  components  of  r  in  the  3  coordinate 
directions  are 


g6r,  a6r,  and  yr 


for  the  x,  y,  and  z  axes  respectively. 

The  error,  e,  or  the  difference  between  the  data  and  this  model  is 

e  «  E  (z±  -  yr.^2  +  (y1  -  a«ri)2  +  (x±  -  66^) 2 

where  i  extends  over  the  time  window  of  interest. 

e  =  Z  x2  +  y2  +  z2  -  2Er^  (yz^  +  aSy^  +  S5x^)  +  Er2 

To  estimate  the  waveforms  at  the  arbitrary  azimuth,  0,  and  incident 

angle,  <f>,  take  the  partial  derivative  of  e  with  respect  to  r  ,  the  n^1  sample 

n 

in  the  time  series  representation  of  the  model. 


—  =  2(yz  +  a5y  +  BSx  )  +  2r 
dr  n  n  n  n 

n 

Setting  this  to  zero 


r  ■  yz  +  a5y  +  05x 
n  n  Jn  n 
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Substituting  the  r  back  into  e, 
n 

e  =*  Zx2  +  y2  +  z2  -  2Ir2  +  Zr2 
=  Ex2  +  y|  +  -  Er2 

To  optimize  r  with  respect  to  0  and  41  we  must  minimize  e.  Thus  we  must 
maximize 


But 


Er2  =  E  (yz  +  otSy  +  B5x)2  =  v 

where  v  is  the  same  function  we  optimized  in  our  earlier  development. 

Thus  our  maneuver  to  avoid  explicitly  invoking  the  straight-line  back- 
and-forth  P-wave  model  gains  us  nothing.  The  process  of  finding  the  azimuth 
and  angle  from  the  vertical  at  which  the  mean  square  excursion  is  maximized 
is  identical  to  the  process  of  estimating  the  azimuth  and  incident  angle  of 
a  P-wave  model  of  straight-line,  back-and-forth  motion. 
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A  Surface-Wave  Processor 


The  surface-wave  processor  described  in  the  Teledyne  Geotech  report 
SDAC-TR-77-14  has  been  extended  in  this  study;  the  ellipticity  ratio  is  no 
longer  a  required  input  parameter  but  is  estimated  simultaneously  with  the 
azimuth  and  spectra.  This  feature  has  proved  necessary  since,  in  practice, 
we  have  not  known  the  ellipticity  ratios  for  the  stations  and  frequencies 
(and  azimuths)  of  the  data  we  have  processed  and  have  had  to  guess  them. 

The  development  of  the  single  mode  surface-wave  processor  with  unre¬ 
strained  ellipticity  follows,  see  Figure  3.  (We  shall  see  that  a  simple 
processor  which  uses  only  horizontal  components  and  maximizes  the  transverse 
to  residual  energy  ratio  gives  better  results.) 

We  may  represent  the  anticipated  fundamental  mode  surface-wave  by  its 
three-component  particle  motion: 

X  is  the  complex  Fourier  coefficient, 

®  jjj/t 

at  the  n  frequency,  of  the  vertical 
component  of  Rayleigh  wave 
displacement 

ie  A.  is  the  corresponding  coefficient 

n  n 

of  the  radial  component,  which  leads 
the  vertical  by  90  degrees  of  phase 
angle  and  differs  from  it  in 
amplitude  by  a  factor  (the 
ellipticity  ratio) 

L  is  the  coefficient  of  the  Love  wave 

n 

displacement,  transverse  to  the 
Rayleigh  displacement  and  unrelated 
to  it  in  both  phase  and  amplitude. 


If  the  azimuth  along  which  this  signal  model  is  propagating  is  *,  then 
the  three-components  we  expect  to  record  are 


vertical: 

north: 

east: 

If  we  let 


K 


■it  4  cos*  -  £  sin* 
n  n  n 


■it  4  sin*  +  £  cos* 
n  n  n 


a  =  cos*  3  =  sin* 

then  the  sum  of  the  squares  of  the  difference  between  the  data  and  the 
signal  model  may  be  written 


E<4  ,  £, 
m  m 


»  Z 
n 


-  4 


l2+  I: 


-  (-te  4  a  -  £  3)  1 2 
n  n  n  ' 


-  (-c  t  4  3 
n  n 


+  £  a) 
n 


The  F-statistic  for  the  model  is  given  by 


F<V  V  ♦>  *  2 


h  <* +  OKI2  +  i«„i2 


«V  V  *> 


We  may  minimize  E  with  respect  to  4  and  £  ,  the  model  spectra.  In 

n  n 

the  usual  fashion  by  taking  the  partial  derivatives  of  E  with  respect  to 
each  of  the  real  and  imaginary  parts  of  each  of  the  spectral  coeff icients. 
Setting  each  derivative  equal  to  zero  we  get 


4  *  (1  +  e2)  *  [z  -  -Ce  (ay  +  Sx  )] 

n  n  n  n  n  n 


£  »  ax  -  3y 
n 
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Substituting  these  results  into  E,  the  sum  of  the  squares  of  the  differ¬ 
ence  between  the  signal  model  and  the  data,  we  get 


E (<fc)  =  El ■Lc  z 
n1  n  n 


(ayn  +  (5xn)|2  (1  +  c2)"1 


If  we  assume  now  that  e  is  constant  over  the  range  of  the  summation, 

n 

that  is, 

e  =  e,  a  constant  with  respect  to  n 
n 

we  may  express  E  as 

(1  +  e2)E  =  e2h  +  a2g  +  S2d  -  2agc  -  2e(aa  +  fSb) 

in  which 


a  =  -?  I  (y*z  -  y  z*) 
2  n  ■'n  n  n  n 


b  2  j  I  (x*z  -  x  z*) 
2  n  n  n  n  n 


c  =  -  ts-  E  (y*x  +  x*y  ) 
2  n  n  n  n  n 


d  =  £ | x  |2 
n1  n1 


8  5  SlyJ2 


h  =  I 

n 


z 

n 


2 
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We  now  take  the  partial  derivatives  of  E  in  this  form,  first  with 
respect  to  e,  and  then  with  respect  to  $,  and  set  each  result  to  zero. 
Taking  the  partial  derivative  of  E  with  respect  to  e  and  setting  it  to 
zero  yields 

e  E  =  eh  -  (aa  +  0b) 

Substituting  this  expression  back  into  the  previous  equation,  we  eliminate 
E  with  the  result  that 

A  e2  +  Be  -  A  =  0 

in  which 

A  =  aa  +  0b 

B  =  h  -  (a2g  -  2a0c  +  02d) 

Now,  taking  the  partial  derivative  of  E  with  respect  to  $  and  setting  it  to 
zero  yields 

eU  =  -V 


in  which 


U  =  ab  -  0a 

V  =  a0(g  -  d)  +  (a2  -  02)c 
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Substituting  this  expression  into 


Ae2  +  Be  -  A  =  0 


yields 


V (AV  -  BU)  -  AU2  =  0 

Finally,  substituting  for  A,  B,  U,  and  V,  we  have  that 

a3W  +  aB2X  +  ag2Y  +  S3Z  =  0 


in  which 


W  =  a(c2  -  b2)  +  (g  -  h)bc 

X  =  ac(g  +  h  -  2d)  +  b(g  -  h)(g  -  d)  -  b(b2  +  c2  -  2a2) 

Y  =  bc(d  +  h  -  2g)  +  a(d  -  h) (d  -  g)  -  a(a2  +  c2  -  2b2) 

Z  =  b(c2  -  a2)  +  (d  -  h)ac 
This  may  be  written  as 

Z  tan3<j>  +  Y  tan2<t>  +  X  tan<()  +  W  =  0 

which  is  a  cubic  equation  in  tan<f>.  The  cubic  equation  is,  of  course,  solved 
analytically  without  computer  search.  This  completes  the  development. 
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As  noted  above,  the  purpose  of  this  study  was  to  develop  a  regional 
event  processor  based  upon  the  fundamental-mode  surface-wave  technique 
developed  by  Smart  (1977).  The  study  was  successful  and,  among  other  things, 
it  has  enabled  us  to  explain  the  past  success  of  the  technique  as  an 
azimuth  estimator  when  applied  to  phases. 

Test  Data  Set 

Thirty-one  events  recorded  at  RK.ON,  well  distributed  azimuthally  about 
the  station  and  covering  the  spectrum  of  signal-to-noise  ratios,  were 
selected  as  a  test  data  set  for  this  study  (Figure  4  and  Table  I).  Signal 
and  noise  spectra  have  been  computed  for  all  the  phases.  Figure  5  shows 
the  spectra  of  four  events,  the  two  highest  frequency  and  the  two  lowest 
frequency  spectral  maxima  in  the  data  set.  (In  general,  the  low-frequency 
signals  were  those  which  traversed  complex  mountainous  regions  on  their  way 
to  RKON,  e.g.,  events  from  California,  Oregon,  Washington,  and  British 
Columbia. ) 

Figure  5  also  shows  composite  noise  spectra  by  seasons.  Contrary  to 
expectations,  the  noise  spectra  were  found  not  to  vary  significantly  over 
the  seasons.  As  can  be  seen  in  the  figure,  the  noise  peaks  below  0.4  Hz, 
whereas  all  signal  spectra  peak  at  or  above  that  frequency.  This  obser¬ 
vation  determined  the  lower  edge  of  the  frequency  window  chosen  for  L 

O 

processing. 

The  Analysis 

For  the  analysis  of  waves  the  modified  version  of  the  Smart 
processor,  which  estimates  particle  motion  ellipticity  as  well  as  azimuth 
and  waveform,  was  used.  Previously,  ellipticity  has  been  an  input  parameter. 
The  frequency  band  examined  was  0.5  and  3.0  Hz.  The  windows,  51.2  seconds 
long,  were  hand  picked  by  an  analyst.  Good  back  azimuths  were  obtained  in 
general:  the  average  of  the  absolute  errors  was  6.7  degrees  (11.1  degrees 
rms).  Some  of  the  low-frequency  events  yielded  the  least  accurate  estimates. 

A  variety  of  parameters  estimated  by  the  L  process  was  monitored 
during  the  experiment  in  a  search  for  a  suitable  detection  index,  i.e.,  one 
which  separates  signals  from  noise.  The  parameters  investigated  were: 
the  F-statistic,  the  mean  recurrence  period,  and  the  derivative  92/E3<f>2, 


Figure  4.  Epicenters  of  31  regional  seismic  events  recorded  at  Red  Lake, 
Ontario  (RKON).  These  3-component  records  form  the  test  data 
set  for  this  report  (sec  Table  I). 
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REGIONAL  EVENTS  AT  RKON 
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explosion 


LOG  POWER 


Figure  5.  (a)  Composite  noise  spectra  for  RKON,  for  four  seasons 

(b)  Peak-normalized  spectra  of  the  lowest-  and  highest- 
frequency  event  nos.  1,  3,  13  and  16  (see  Table  I) 
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where  E  =  error  and  <p  =  back  azimuth  (see  Smart,  1977).  The  rms  excursions 
for  the  model  in  the  vertical,  radial,  and  transverse  directions,  denoted 
by  Z,  R,  and  T,  were  also  examined. 

The  derivative  did  not  follow  a  discernible  pattern.  The  F-statistic 
and  mean  recurrence  period  were  consistently  low  for  both  noise  and  signal 
and  did  not  separate  them  (see  Figure  6  for  F-statistic  histogram).  The 
ratio  T/R  proved  a  somewhat  useful  discriminator,  but  its  separation  of 
the  populations  was  not  complete  enough  to  qualify  it  as  a  detector 
(Figure  7).  It  was  observed  that  T/R  <  2.0,  when  it  occurred,  accompanied 
poor  azimuthal  estimates,  i.e.,  error  >  10  degrees,  and  that  the  smaller 
T/R,  the  larger  the  azimuthal  error. 

These  observations  suggested  that,  by  itself,  the  figure  traced  by 

the  particle  motion  in  the  horizontal  plane  might  serve  to  indicate  surface- 

wave  azimuth.  An  algorithm  was  coded  to  find  in  each  time  window  that 

orientation  in  the  horizontal  plane  where  the  rms  excursion,  in  the 

frequency  band  of  interest,  was  minimized.  The  azimuth  of  the  orientation 

was,  on  the  average,  about  as  good  an  estimate  of  surface  wave  origin  as  the 

back  azimuth  from  the  original  coherent  processor.  Moreover,  it  was  more 

stable.  In  low  signal-to-noise  records,  where  T/R  fell  below  2.0,  the 

estimated  azimuth  was  more  accurate  than  that  of  the  coherent  processor. 

All  this  appears  to  explain  why  the  Smart  processor  has  been  successful  for 

estimating  L  azimuths,  but  unsuccessful  for  detecting  L  arrivals.  The 
8  8 
azimuth  estimations  have  been  controlled  by  the  dominant  transverse  Love 

component,  but  the  F-statistics  have  been  kept  low  by  the  incoherence 

between  the  vertical  and  radial  components.  This  incoherence  results  from 

the  mixing  of  fundamental  Rayleigh  modes  with  higher  ones. 

The  results  suggested  another  use  for  this  same  simple  algorithm 
which  monitors  the  elongation  of  the  figure  described  by  particle  motion, 
that  is,  to  flag  the  arrival  of  P-waves  and  indicate  their  azimuths.  6.4 
second  time  windows — an  order  of  magnitude  smaller  than  those  for  the  L 
waves — were  used  to  search  in  the  vicinity  of  the  expected  P  arrivals. 

The  principal  axis  of  the  elongated  figure  was  taken  to  indicate  back 
azimuth,  and  the  amount  of  elongation,  that  is,  the  ratio  of  the  principal 
long  axis  to  the  short  axis,  was  taken  as  the  index  of  signal  presence. 

Taking  that  ratio  at  2.4  as  the  detection  threshold,  in  the  band  0.9  to  6.0  Hz, 
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FREQUENCY  OF  OCCURRENCE 


Figure  6.  F-statistics  for  the  noise  and  for  the  surface  waves 
of  the  present  test  data  set  (see  Table  I)  for  the 
frequency  band  of  0.5  to  3.0  Hz  show  negligible  separation 
of  populations. 


80 %  of  the  signals  were  detected  with  a  high  false  alarm  rate  of  28  per  hour. 
The  average  azimuthal  error  was  6.0  degrees,  and,  in  the  case  of  the  low- 
frequency,  low  signal-to-noise  events,  these  P  azimuthal  estimates  were 
more  accurate  than  those  from  the  L  waves.  Thus,  in  spite  of  the  high 
L^/P  amplitude  ratios,  the  P  wave  appeared  as  useful  for  signal  detection 
and  for  azimuth  estimation  as  L^,  within  the  restraints  of  the  present 
investigation. 

However,  this  technique  made  use  of  only  two  components  of  the  data, 

and  in  the  next  set  of  tests  the  three-component  P-wave  process,  developed 

above,  proved  more  useful.  Since  it  estimates  the  emergence  angle  of  the 

wave  and  thus  takes  account  of  the  motion  up  and  away  from  the  source,  P-wave 

particle  motion  in  three-components  is  free  of  the  180  degree  ambiguity  of 

azimuth  inherent  in  L  particle  motion  (and  in  the  observations  of  P  in 

8 

the  horizontal  plane  alone,  as  above) .  The  average  error  in  estimated 
azimuth  (from  the  true  geodesic  azimuth)  over  the  entire  data  set  is  7.0 
degrees  for  the  3-component  P  processor  (9.1  degrees  rms) ,  Setting  an  F 
threshold  which  rejects  only  13  percent  of  the  estimates  reduces  the  average 
azimuth  error  to  less  than  6.0  degrees  (7.4  degrees  rms). 

Since  P  waves  may  be  emergent  as  well  as  impulsive  (see  Figures  8  and 
9),  the  three-component  P  process  was  not  simply  applied  to  the  one  arrival 
window  chosen  by  the  analyst,  but  was  allowed  to  search  that  vicinity 
with  a  sliding  window.  Of  the  several  resulting  estimates,  those  with 
unacceptably  low  emergence  angles  (less  than  20  degrees  from  the  horizontal) 
were  rejected  and  the  azimuth  was  picked  from  those  remaining  on  the  basis 
of  the  F-statistic. 

Combining  the  estimated  L  azimuth  with  that  of  the  P  wave,  simply 

8 

by  taking  their  mean,  increased  the  accuracy.  The  average  difference 
between  the  mean  estimated  azimuth  and  the  true  geodesic  azimuth  was  4.9 
degrees  (7.0  degrees  rms).  Moreover,  the  F-statistic  computed  by  these 
processors  serves  to  separate  poor  azimuthal  estimates  from  the  population: 
the  azimuthal  estimates  (over  80%  for  this  data  set)  which  passed  the 
arbitrary  F  threshold  set  for  this  study  differed  from  the  known  geodesic 
azimuths  by  an  average  of  3.9  degrees  (4.9  degrees  rms).  It  is  the  author's 
experience  that  such  accuracy  is  comparable  to  or  better  than  that  from 
a  well-designed  array  of  sensors. 
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Figure  9. 
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EVENT  17 


An  emergent  P-wave  signal  selected  from  the  test  data 
set  of  this  study  (see  Table  I,  event  17). 
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Discussion 


The  limitations  of  this  study  have  not  permitted  the  conclusion  of 

the  research.  The  optimum  frequency  band,  the  optimum  signal  flag,  and 

the  most  useful  combination  of  P  and  L  information  have  not  been 

g 

determined.  Two  additional  detection  criteria  have  suggested  themselves 
during  this  work.  One  is  the  observed  stability  of  orientation,  over 
several  time  windows,  of  the  elongated  particle-motion  envelope,  both 
in  the  P-wave  portion  of  the  record  and  in  the  portion.  The  other  is 
the  observed  ninety  degree  rotation  of  the  estimated  azimuth  as  the  P-wave 
passes  and  the  wavetrain  arrives.  These  points  demand  further  study. 

The  parameters  estimated  by  the  processors  have  not  been  fully 

exploited  as  yet.  The  emergency  angle  for  P  waves  has  been  used  here,  and 

to  good  advantage,  to  avoid  confusing  waves  with  P  waves.  They  both 

have  back-and-forth  straight-line  motion,  but  P  motion  is  well  out  of  the 

horizontal  plane,  while  L  motion  is  largely  confined  to  it.  Thus,  the 

g 

estimated  emergency  angle  is  an  important  criterion  for  P  processing. 

But  the  ellipticity  measurements  from  the  surface-wave  processor  have 
yet  to  be  employed.  In  this  study,  of  RJKON  data,  ellipticity  clusters 
around  0.6,  and  the  outliers  are  associated  with  inaccurate  azimuths. 

Now  that  it  is  known  to  lie  around  0.6  for  RKON,  the  ellipticity  can  be 
held  fixed,  which  will  improve  detection  sensitivity,  and  may  also  bring 
those  outliers  in  closer. 

There  are  grounds  for  suspecting  that  the  long  time  windows  employed 
for  processing  may  have  had  an  adverse  effect  upon  the  results  of 
this  study,  at  least  as  far  as  detection  is  concerned.  Moreover,  the 
three-component  P  processor  has  not  been  applied  to  noise  to  evaluate 
its  sensitivity  as  a  detector.  So,  further  trials  using  various  window 
lengths  are  in  order. 

Finally,  it  appears  from  the  research  carried  out  thus  far,  that 
the  principal  utility  of  particle-motion  processing  lies  not  in  signal 
detection,  where  the  simple  power  detector  is  successful,  but  in  azimuth 
and  distance  determination.  Of  the  missed  signals  in  the  P-wave  detection 
trials  discussed  earlier,  the  P-wave  processor  estimated  the  back  azimuth 
of  half  of  them  to  within  five  degrees  of  their  true  back  azimuth,  even 
though  it  did  not  "detect"  them.  One  of  these  signals  was  not  visible  even 


-34- 


to  the  analyst.  The  time  between  P  and  L^,  measured  from  the  power  rise 
at  the  signal  onset  down  to  the  point  where  the  particle-motion  orientation 
rotates  90  degrees,  yields  a  distance  estimate.  In  the  case  of  P-waves 
from  sources  at  distances  greater  than  20  degrees  from  the  recording 
station,  the  emergency  angle  also  provides  an  estimate  of  source  distance. 
Particle-motion  processing  can  also  be  exploited  for  picking  phases. 

All  of  this  potential  must  be  addressed  in  future  research. 
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